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IC 4663: The first unambiguous [WN] Wolf-Rayet central 
star of a planetary nebula* 
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ABSTRACT 

We report on the serendipitous discovery of the first central star of a planetary nebula 
(PN) that mimics the helium- and nitrogen-rich WN sequence of massive Wolf-Rayet 
(WR) stars. The central star of IC 4663 (PN G346.2-08.2) is dominated by broad He II 
and N V emission lines which correspond to a [WN3] spectral type. Unlike previous 
[WN] candidates, the surrounding nebula is unambiguously a PN. At an assumed 
distance of 3.5 kpc, corresponding to a stellar luminosity of 4000 Lq, the V — 16.9 mag 
central star remains 4-6 mag fainter than the average luminosity of massive WN3 stars 
even out to an improbable d — S kpc. The nebula is typical of PNe with an elliptical 
morphology, a newly discovered Asymptotic Giant Branch (AGB) halo, a relatively 
low expansion velocity (wexp = 30 km s""'^) and a highly ionised spectrum with an 
approximately Solar chemical abundance pattern. The [WN3] star is hot enough to 
show Ne VII emission (T* = 140±20 kK) and exhibits a fast wind (woo = 1900 km s^^), 
which at d = 3.5 kpc would yield a clumped mass loss rate of M = 1.8 x 10~*Mo yr~^ 
with a small stellar radius (i?* =0.11 Rq). Its atmosphere consists of helium (95%), 
hydrogen (< 2%), nitrogen (0.8%), neon (0.2%) and oxygen (0.05%) by mass. Such an 
unusual helium-dominated composition cannot be produced by any extant scenario 
used to explain the H-deficiency of post-AGB stars. The 0(He) central stars share 
a similar composition and the discovery of IC 4663 provides the first evidence for 
a second He-rich/H-deficient post-AGB evolutionary sequence [WN]^'0(IIe). This 
suggests there is an alternative mechanism responsible for producing the majority 
of H-deficient post-AGB stars that may possibly be expanded to include other He- 
rich/H-deficient stars such as R Coronae Borealis stars and AM Canum Venaticorum 
stars. The origin of the unusual composition of [WN] and 0(He) central stars remains 
unexplained. 

Key words: planetary nebulae: individual: PN G346.2— 08.2 - planetary nebulae: 
general - stars: wolf-rayet - stars: AGB and post-AGB - stars: abundances - stars: 
mass-loss 



1 INTRODUCTION 

Planetary nebulae (PNe) are the result of low-intermediate 
* Based on observations made with Gemini South under pro- m^ss stars (-1-8 Mq) that have experienced heavy mass 

gramme GS-2011A-Q-65 (PI: B. Miszalski). loss during the asymptotic giant branch (AGB) phase. 

t E-maih brent@saao.ac.za Their gaseous envelopes are ionised by hot pre-white dwarfs 
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(Teff~30-140 kK) whose diverse spectroscopic appearance 
(Mendez 1991) reflects different wind properties and at- 
mospheric compositions. A subset of central stars have H- 
deflcient atmospheres composed primarily of helium, carbon 
and oxygen with fast and dense winds that produce broad 
emission lines. Their spectroscopic appearance mimics the 
carbon-sequence of massive Wolf-Rayet (WR) stars and they 
are classified into either [WO] or [WC] types, the brack- 
ets distinguishing them apart from massive stars (van der 
Hucht 2001), depending on which emission lines are present 
and their relative intensities (Crowther, De Marco & Bar- 
low 1998; Acker & Neiner 2003). Despite recent surveys that 
have dramatically increased the number of central stars that 
belong to the [WC] sequence (Acker & Neiner 2003; Gorny 
et al. 2004, 2009, 2010; DePew et al. 2011), their origin and 
evolution remain poorly understood. 

Werner & Herwig (2006) reviewed possible scenarios 
that may explain the H-deficiency of [WC] central star^ and 
their progeny the PG1159 stars (Wesemael, Green & Liebert 
1985). After the central star has left the AGB phase a late 
thermal pulse (LTP) or very-late thermal pulse (VLTP) 
would reignite helium-shell burning (e.g. Blocker 2001; Her- 
wig 2001). The star is sent back to the AGB phase, i.e. it 
becomes "born-again" (Schonberner 1979), and any residual 
hydrogen is either burned or diluted below the observable 
detection limit. The AGB final thermal pulse (AFTP) vari- 
ant may also be invoked to explain some instances where a 
small amount of observable hydrogen remains. At this point 
the star enters the [WC]-sequence with an H-deficient at- 
mosphere mostly made of helium (~30-50%) , carbon (~30- 
60%) and oxygen (~2-20%). 

A post-AGB H-deficient evolutionary sequence 
[WCL]-^[WCE]->PG1159 has been proposed in which 
late-type [WC] central stars evolve into early-type [WC] 
and [WO] central stars, followed by PG1159 stars which 
appear preferentially in the oldest PNe (Napiwotzki & 
Schonberner 1995) and in isolation once the surrounding 
PNe have dissipated into the interstellar mediumlj Unless 
there is a residual amount of hydrogen the end product will 
be a non-DA white dwarf. This evolutionary sequence is 
supported by the strong similarities in derived atmospheric 
compositions (Werner & Herwig 2006) and some transition 
objects (e.g. Bond 2008). A key question is whether the H- 
deflciency of all members of the sequence can be explained 
by the LTP, VLTP and AFTP scenarios. Corny & Tylenda 
(2000) found the evolutionary status of nebulae surrounding 
H-deflcient central stars to be consistent with H-rich central 
stars, rather than following the paths predicted by LTP 
or VLTP scenarios. It is therefore debatable whether the 
few documented LTP and VLTP events (Schonberner 
2008) mandate the same origin for the larger population 
of H-deflcient central stars. Indeed, there may even be an 
alternative common-envelope (De Marco & Soker 2002; De 



^ Massive WC stars are extremely H-deficient, but for a differ- 
ent reason. A strong wind has peeled off the outer H-containing 
layers exposing in turn the C- and O-rich products of triple-o 
He-burning. 

^ This differs from massive WC stars which appear to have their 
subtype determined by their initial metallicity, although some 
evolution from late to early may also occur (Smith & Maeder 
1991). 



Marco 2008; Hajduk, Zijlstra & Gesicki 2010) or nova-like 
explanation for such objects (Wesson, Liu & Barlow 2003; 
Wesson et al. 2008; Lau, De Marco & Liu 2011). 

Even more difficult to explain are the helium-rich H- 
deficient stars which do not fit the carbon-rich post-AGB 
H-deficient sequence. These include four 0(He) stars (Ranch 
et al. 1998), two of which are central stars of PNe (Ranch 
et al. 1994, 1996), R Coronae Borealis (RGB) stars (Clay- 
ton 1996), extreme helium B stars and helium-rich subdwarf 
O stars (see e.g. Werner & Herwig 2006). With typical he- 
lium mass fractions >90% and trace amounts of other ele- 
ments the chemical compositions of these stars have never 
been predicted by the born-again scenarios, i.e. they sim- 
ply should not exist! On the whole these stars are scarce 
compared to the carbon-rich sequence and probably form 
a second, parallel post-AGB H-deficient sequence. Evidence 
for such a sequence is tentative however, based primarily on 
their common composition rather than other evolutionary 
ties. It has been suggested that RGB stars are the progeni- 
tors of 0(He) stars, although the circumstances of this pos- 
sible connection remain unclear. After the discovery of a low 
i6q^i8q j.g^^-Q ^jy Clayton et al. (2007) strongly pointing to a 
merger origin for RGB stars, an end to the protracted debate 
on the formation of RGB stars seemed to be in sight. How- 
ever, not only is the RGB evolutionary path complicated by 
a connection with [WC] central stars (Clayton & De Marco 
1997; Clayton et al. 2006; Clayton et al. 2011a), but there 
are indications of a mixed origin within the RGB class itself 
(Clayton et al. 2011b). 

It is imperative that more members of the He-rich/H- 
deficient class are found and studied to establish a second 
evolutionary sequence. If this can be achieved, then there 
should exist another mechanism responsible for H-deficiency, 
one which may also help explain the origin of the carbon- 
rich sequence. A long sought after solution may be [WN] 
type WR central stars that mimic the helium- and nitrogen- 
rich WN sequence in massive stars (Smith, Shara & Moffat 
1996) . Although some candidates have been proposed, none 
have turned out to be satisfactory and the existence of the 
class remains unproven. 

Here we report on the serendipitous discovery of a [WN] 
central star in the PN IG 4663 (PN G346.2-08.2, Fleming 
& Pickering 1910). Observations of IG 4663 were made dur- 
ing the Gemini South programme GS-2011A-Q-65 (PI: B. 
Miszalski) which aimed to discover new close binary central 
stars of PNe via radial velocity variability (e.g. Miszalski 
et al. 2011a). Program targets were selected based on the 
nebula trends identified by Miszalski et al. (2009b) to be 
provisionally associated with binarity that includes bipolar 
nebulae, nebulae with coUimated outfiows (jets) and low- 
ionisation structures (e.g. Gongalves et al. 2001). As noted 
by Miszalski et al. (2009b), a large number of PNe with WR 
central stars overlap with this morphological prescription 
and so it was not unexpected that new examples would be 
found during the programme. 

This paper is structured as follows. Section [2] reviews 
previous [WN] candidates. Section|3]describes Gemini South 
spectroscopy and atmospheric modelling of the previously 
unobserved central star of IG 4663. The basic nebula prop- 
erties and chemical abundances are presented in Section 2] 
Section [5] establishes the [WN] nature of the central star by 
proving the bona-fide PN nature of IG 4663 and introduces 
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the second post-AGB H-deficient evolutionary sequence. We 
conclude in Section [6] 



2 PREVIOUS [WN] CANDIDATES 

Around 8% of massive Galactic WN stars are known to 
have surrounding nebulae (Stock & Barlow 2010). Until im- 
proved observations could reliably determine the luminosity 
of their central stars, some of these nebulae were misclassi- 
fied as PNe. As the distances to PNe are extremely difficult 
to estimate this has not helped clarify the status of many 
candidates. Figure [T] shows M 1-67 which was included in 
the PN catalogue of Perek & Kohoutek (1967), however its 
WN8h star WR124 was only later proven to be a massive 
star by Crawford & Barlow (1991a) based on the kinematic 
distance measured from the interstellar Na 1 D2 absorption 
line|j The flocculent appearance of the nebula, evident in 
Fig. [ij but best seen in HST observations (Grosdidier et al. 
2001a), also bears more resemblance to nebulae around mas- 
sive WN stars including RCW 58 around WR40 (Chu 1982) 
and NGC6888 around WR136 (Wendker et al. 1975). An- 
other similar case to M 1-67 is We 21 (Duerbeck & Reipurth 
1990) which was proven to host a massive star by Crawford 
& Barlow (1991b). 

Other [WN] candidates have proven very difficult to es- 
tablish as promising candidates for clear-cut low-mass ana- 
logues of the massive WN sequence. The most extensively 
studied is LMC-N66 in the Large Magellanic Cloud (see e.g. 
Pefia et al. 1994, 1997, 2004, 2008; Hamann et al. 2003). Mul- 
tiple outbursts have been recorded in the central star which 
shows a mid-WN spectral type during outburst. Pefia et al. 
(2004) studied the surrounding nebula (Fig. [l]) which resem- 
bles point-symmetric PNe in the Milky Way (e.g. Fleming 
1, Lopez et al. 1993) and shows multiple emission compo- 
nents that span 130-160 km s~^. The luminosity spans log 
L/Lq =4.6-5.4 (Hamann et al. 2003), too high to be con- 
sistent with a typical PN central star. Hamann et al. (2003) 
presented a number of possible explanatory scenarios, none 
of which are without difficulty, although a massive star in 
a common-envelope phase with a low mass companion or a 
massive white dwarf (WD) accreting matter in a close bi- 
nary seem to be preferred. Binary stellar evolution is there- 
fore likely responsible for the WN features of LMC-N66, as 
is similarly seen in symbiotic novae such as RR Tel and PU 
Vul (Thackeray & Webster 1974; Nussbaumer 1996). 

Morgan, Parker & Cohen (2003) discovered a WN6 star 
in the nebula PM5 and argued for a PN status based on 
the nebula density, size, ionised mass, morphology, near- 
infrared and mid-infrared properties mostly being similar 
to PNe. A non-PN status is supported by the anomalously 
high expansion velocity (165 km s~^), low Galactic latitude 
(fo = +0.69 deg), high reddening {E{B -V) = 3.0-3.5 mag) 
and the analysis by Todt et al. (2010b). Furthermore, Fig.[T] 
shows a previously unpublished VLT FORS2 (Appenzeller et 
al. 1998) image observed under ESO programme ID 077. B- 
0430(B) (PL M. Pefia). The 20 sec image taken through the 
H_alpha-|-83 filter (central wavelength (Accn) and full-width 



^ Marchenko, Moffat & Crowther (2010) recently confirmed this 
result by analysing the nebula expansion parallax with 2-epoch 
HST observations. 



at half-maximum (FWHM) of 656.3 nm and 6.1 nm, respec- 
tively) reveals a flocculent appearance which was not visible 
in the lower resolution SuperCOSMOS Ha Survey discovery 
images (SHS, Parker et al. 2005). Similarly, the WN6 central 
star of A 48 (AbeU 1966; Wachter et al. 2010; DePew et al. 
2011), is also likely to be a massive star, however this object 
has not yet been studied in detail. It is also located at a low 
Galactic latitude {b = -f-0.4 deg), has a double- ring mor- 
phology (Stock & Barlow 2010) and is also quite reddened 
{E{B -V) = 2.0 mag, DePew et al. 2011). 

Several other candidates have also been identified based 
on the presence of weak nitrogen emission features in addi- 
tion to their [WO/WC] spectra. The [W04/WC4] central 
stars of NGC 6751 (AUer 1968; Koesterke & Hamann 1997), 
PMRl and PM3 (Morgan, Parker & Russeil 2001; Parker 
& Morgan 2003) exhibit weak N fV and N V emission fea- 
tures. Further examples are evident amongst the [W04pec] 
and [WO] central stars studied by Acker & Neiner (2003). 
In the case of NGC 6751 nitrogen accounts for no more than 
1% by mass in an otherwise typical [WO/WC] atmospheric 
composition of 55% helium, 30% carbon and 15% oxygen 
(Koesterke & Hamann 1997). Others were analysed by Todt 
et al. (2008) who found nitrogen contributed no more than 
1-2% by mass. 

Todt et al. (2010a, 2010b) found the central star of PB 8 
to show clear nitrogen emission features and determined a 
nitrogen mass fraction of 2%. Its overall atmospheric compo- 
sition was found to be made of oxygen and carbon each 1.3% 
by mass, helium 55% and hydrogen an unusually high 40%. 
This is most unlike [WO/WC] central stars which have sig- 
nificantly higher fractions of carbon and oxygen (Crowther 
2008). The PN nature of PB 8 is demonstrated by its rel- 
atively low expansion velocity (19 ± 6 km s~^) and espe- 
cially our discovery in Fig. [T] of a faint AGB halo (Corradi 
et al. 2003). Todt et al. (2010a) proposed that PB 8 belongs 
to a [WN/WC] classification, an analogue of the massive 
WN/WC stars (e.g. Massey & Grove 1989; Crowther et al. 
1995a) which account for <5% of massive WR stars (van 
der Hucht 2001). Therefore, strictly speaking, PB 8 is not a 
true [WN] nucleus due to the additional [WC] features. 

To summarise, no [WN] candidate has been unambigu- 
ously proven to be surrounded by a PN. While PB 8 is cer- 
tainly a PN, it appears to be a [WN/WC] hybrid rather than 
a strict [WN]. 



3 THE CENTRAL STAR 

3.1 Observations 

We used the Gemini Muhi-Object Spectrograph (CMOS, 
Hook et al. 2004) on Gemini South to obtain three spec- 
tra of the central star of IC 4663 on 27 June, 2 July and 
3 July 2011 as part of the GS-2011A-Q-65 programme (PI: 
B. Miszalski). The B1200 grating was used with a 0.75" slit 
oriented at a position angle of 90° to provide wavelength 
coverage from 3767-5228 A (27 June and 3 July) and 4507- 
5977 A (2 July) at 1.6 A resolution (FWHM) with a recip- 
rocal dispersion of 0.23 A pixel"^ (1x2 binning in disper- 
sion and spatial axes, respectively). The science exposure 
times were 1800 s (27 June) and 2400 s (2 and 3 July) and 
were accompanied by contemporaneous arc lamp exposures. 
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Figure 1. A selection of previous [WN] candidates. The M 1- 
67 (Drew et al. 2005; Gonzalez-Solares et al. 2008) and PM 5 
(see text) images were both taken with an Ha+[N II] filter, while 
LMC-N66 was observed with the clear filter of HST STIS. The 
colour-composite image of PB 8 is made from SHS Ha+[N II] 
(red), SHS Short-Red (green) and SuperCOSMOS Sky Survey 
Bj (blue, Hambly et al. 2001). It reveals for the first time a faint 
AGB halo with a diameter up to ~80" surrounding PB 8 (the 
bright central source). Image dimensions are 120 X 120 arcsec^ 
(M 1-67), 6x6 arcsec^ (LMC-N66), 40 x 40 arcsec^ (PM 5) and 
4x4 arcmin^ (PB 8). In all cases North is up and East to left. 



The data were reduced with the Gemini iraf package af- 
ter which nebula-subtracted one-dimensional spectra were 
extracted using the APALL iraf task. A non-uniform spatial 
profile in the nebula emission lines produced a small amount 
of over-subtraction at the Balmer lines, He II AA4686, 5412, 
and He I A5876. The one-dimensional spectra were flux cali- 
brated with a spectrum of the spectrophotometric standard 
star EG274 in the usual fashion. 

The V magnitude of the central star is not well deter- 
mined in the literature with discordant estimates of V = 
15.161^23 mag (Shaw & Kaler 1989) and B = 18.0±(0.5- 
1.0) mag (Tylenda et al. 1991). We retrieved the processed 
and combined HST F555VK image (approx. Johnson V 
of total duration 1020 s from the Hubble Legacy Archiv( 
(HLA) and performed aperture photometry on the central 
star following the WFPC2 Photometry Cookbook|j Assum- 
ing V — F555VK = 0.05 mag, suitable for an 05V star, we 
found V — 16.90 ± 0.05 mag and fixed the absolute fiux 
scale of our GMOS spectrum to match the corresponding 
Vega flux at A5480 A. 



uum (except for A ^ 4500 A where S/N > 30) and rectified 
for comparison with the NLTE model described in Sect. 13.31 
The spectrum is dominated by He II and N V emission lines 
typical of early-type WN stars (Crowther et al. 1995b) and 
shares the triangular line profiles of Galactic early-type WN 
stars (Smith et al. 1996). The Ne VII emission lines (Werner, 
Ranch & Kruk 2007) were never before seen in massive WN 
stars, whereas the weak O V and O VI emission are uncom- 
mon (e.g. the WN3pec star WR46 shows relatively strong 
O VI A3811-34, Crowther et al. 1995b). The high S/N spec- 
trum rules out the presence of He I and C IV emission lines 
which, combined with the already mentioned features, cor- 
responds to a clear WN3 subtype (Smith et al. 1996). The 
WN3 subtype is not common amongst massive Galactic WN 
stars where only a few examples are known (e.g. Crowther 
et al. 1995b; Hamann et al. 1995; Smith et al. 1996). The 
equivalent width of He II A4686 (~40 A) is comparable to 
massive WN3 stars, whereas its FWHM is slightly lower 
than WN3 stars (~18.5 A). This behaviour is also seen in 
C III A5696 in the comparison between the [WC9] central 
star of BD-l-30 3639 and massive WC9 stars (Crowther et al. 
2006). In the rest of this paper we denote the classification 
as [WN3] to indicate it is a bona-fide central star of a FN 
(later proven in Sect. l5.l"|) rl 

3.3 NLTE modelling 

We have analysed the GMOS spectrum of IC 4663 using 
CMFGEN (Hillier & Miller 1998) which solves the radiative 
transfer equation in the comoving frame subject to sta- 
tistical and radiative equilibrium, assuming an expanding, 
spherically symmetric, homogeneous atmosphere, allowing 
for metal line blanketing and wind clumping. 

Our adopted model atom includes Hi, Hel-ll, Civ, 
Niv-v, Oiv-vi, Neiii-viii, Siiv and Fev-ix. In total 512 su- 
perlevels, 1907 levels, and 28397 non-local thermodynamic 
equilibrium (NLTE) transitions are simultaneously consid- 
ered, with H/He, C/He, N/He, O/He variable, and Si and 
Fe mass fractions fixed at Solar values (Asplund et al. 2009). 
We adopt a Ne mass fraction of 0.2%, based on a number 
ratio of Ne/0 = 0.25 from Wang & Liu (2008) and the Solar 
O mass fraction of 0.55% (Asplund et al. 2009). 

In view of the relatively weak emission line spectrum 
of IC 4663, the supersonic velocity law, «oo(l — R/rf , is 
merged with a subsonic structure calculated using the NLTE 
plane-parallel code TLUSTY (Hubeny & Lanz 1995) Q In or- 
der to reproduce the morphology of Hell A4686 we obtain 
^oo = 1900 kms~^ and /? ~ 2.5, the former being similar 
to iioo = 2000 km s~^ as measured for the massive WN3 
star WR152 (Crowther et al. 1995b). A clumped mass loss 
rate is calculated for an adopted volume filling factor of / = 
0.1, although this is poorly constrained for IC 4663. While 
/3 values close to unity are normally found in weak-wind 



3.2 Spectral features and [WN3] classification 

Figure [2] shows the two 2400 s exposures averaged to achieve 
a signal-to-noise ratio (S/N) of at least ~60 in the contin- 



* http://hla.stsci.edu 

^ http://www.stsci.edu/hst/wfpc2/analysis/wfpc2_cookbook.html 



^ Despite the good match with WN3 stars, we caution that the 
Smith et al. (1996) scheme was solely devised for massive WN 
stars. If a sufficient sample of [WN] stars were to be found, then 
a unified classification scheme may be required (e.g. Crowther et 
al. 1998) which may not necessarily match the existing one. 
^ For our derived T* = 140 kK and logL/L© = 3.6 we adopt 
M = 0.6 Mq resulting in log g = 6.1 for the TLUSTY model 
structure. 
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Figure 2. Rectified Gemini South CMOS spectrum of IC 4663 (black) with the CMFGEN NLTE model (red, see Sect. [331. The main 
emission lines are identified, most belonging to He II and N V, which together with the absence of C IV and He I emission, indicates a 
[WN3] spectral type. The apparent double-peaked profile of He II A4686 is an artefact of nebula subtraction and does not influence the 
NLTE analysis. 



OB stars, larger values implying more extended accelerat- 
ing zones are found in massive WR stars (e.g. Moffat 1996; 
Lepine & Moffat 1999). 

Diagnostic optical lines of He ii A4686 and Ne vii A3890 
enable the stellar temperature of T, — 140±20 kK to be 
obtained. Beyond this range, Nevil weakens at higher tem- 
peratures, where it is replaced by Nevill (e.g. A4340) and at 
lower temperatures, where Nevi strengthens (e.g. A4973). 
High lying recombination lines of Nv (AA4945) and Ovi 
(A5290) serve as reliable N/He and O/He abundance diag- 
nostics, from which mass fractions of 0.8% and 0.05% are 
obtained, respectively. Negligible hydrogen is required to re- 
produce the lower members of the Pickering-Balmer series. 
Hen A5412, Hell-|-H/3 and Hen A4542, permitting an upper 
limit of ~2% by mass. Since there are no suitable diagnostics 
for C/He we investigated whether we could place limits on 
the carbon mass fraction by recalculating the CMFGEN model 
with an elevated carbon abundance. The C IV A5801, 5812 
emission lines became visible for values in excess of 0.1%. 

The CMFGEN model provides a satisfactory fit to the 
CMOS spectrum (Fig. [2|, although Ov A5590 is too weak, 
as is the blend of Nevill A5668 and N v A5670. In addition, 
higher members of the Hen Pickering series are predicted 
to exhibit weak P Cygni profiles, in contrast to observa- 
tions. The theoretical spectral energy distribution (SED) 
was matched to the dereddened CMOS spectrophotome- 
try enabling a determination of luminosity and mass loss 



rate. We obtained an interstellar extinction of E{B — V) 
— 0.35 ± 0.05 (Ay = 1.1 mag) using a standard Galactic 
extinction law (Howarth 1983). This agrees well with the 
value of E{B — V) = 0.4 ± 0.1 mag that Cavichia, Costa & 
Maciel (2010, hereafter CCMIO) measured from the Balmer 
decrement of the surrounding nebula. 

We adopt a distance of 3.5 kpc from the statistical 
Shklovsky distance scale of Stanghellini, Shaw & Villaver 
(2008, hereafter SSV08) which includes Magellanic Cloud 
PNe as calibrators. No individual error was given for the dis- 
tance to IC 4663. We therefore take a conservative approach 
and assume a 2a uncertainty of 60%, which results in a dis- 
tance of 3.5±2.1 kpc. With this distance an absolute magni- 
tude of Afv ~ +3.1 mag is derived for the [WN3] central star 
of IC 4663. The bolometric correction for our CMFGEN model 
is —7.4, resulting in A^boI = —4.3 or log L/Lq — 3.6. For this 
luminosity, a stellar radius of i?, — 0.11 -Rq and a clumped 
mass loss rate of M = 1.8 xlO~*Af0yr~^ are required. As 
indicated in Fig. O the SED of IC 4663 is exceptionally 
hard, with Qo = 10*"='" s~\ Qi = lO*^'"" s-\ Qa = iC^-*'^ 
s~^ for the Hi, He I and Hell continua, respectively. Table[T] 
summarises the stellar parameters of the [WN3] central star 
and includes ranges of distance-dependent quantities for the 
2.1 kpc uncertainty in the distance. 

With the derived stellar properties of IC 4663 we can 
quantitatively judge its wind strength against massive WN3 
stars. Table [2] summarises the stellar properties of massive 
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log Wavelengtli (A) 

Figure 3. Intrinsic and reddened SEDs of the NLTE model (red 
and blue, respectively) with the GMOS spectrum (black). 



Table 1. Stellar parameters of the [WN3] central star of IC 4663. 
Quantities in square brackets represent the range of parameter 
values for the 2.1 kpc uncertainty in distance. 



T. 


140 ± 20 


kK 


logs 


6.1 


(assuming M = 0.6 Mq) 


«oo 


1900 


km s"-"^ 


E{B - V) 


0.35 ±0.05 


mag 


V 


16.90 ±0.05 


mag 


d 


3.5 ±2.1 


kpc (SSV08) 


Mv 


3.1 


mag [2.1-5.1] 


log L/Lq 


3.6 


[2.8-4.0] 


R* 


0.11 


Rq [0.04-0.17] 


M 


1.8 


lO-*M0yr-i [0.5-3.5] 


H 


<2 


mass fraction (%) 


He 


95 


mass fraction (%) 


C 


<0.1 


assumed mass fraction (%) 


N 


0.8 


mass fraction (%) 


O 


0.05 


mass fraction (%) 


Ne 


0.2 


mass fraction (%) 



WN3 stars (Hamann et al. 2006) and IC 4663, together with 
the transformed radius Rt (see Eqn. 1 of Hamann et al. 
2006). The magnitude of Rt diagnoses the wind density (i.e. 
the emission-line strength) which is essentially the same for 
massive WR and [WR] stars of a given spectral type, despite 
the order of magnitude difference in radii (e.g. Crowther et 
al. 2006). Low values of Rt correspond to strong, thick winds 
and high values to weak, thin winds. Since the Rt of IC 4663 
lies between WR152 and WR3, the stellar wind is similarly 
dense to massive WN3 stars, entirely consistent with our 
morphological [WN3] classification. This argues against any 
weak wind classification of the central star of IC 4663 (e.g. 
Of or Of-WR, see Mendez, Herrero & Manchado (1990) and 
Mendez (1991)). 



Table 2. Stellar properties and transformed radii Rt for three 
massive WN3 stars (Hamann et al. 2006) and the [WN3] central 
star of IC 4663. The D quantity is the virind clumping factor. 



Object R*/R0 log L/L(. 



log M D Rt/Re 



WR3 


2.65 


5.6 


2700 


-5.3 


4 


15.1 


WR46 


2.11 


5.8 


2300 


-5.1 


4 


6.8 


WR152 


2.20 


5.3 


2000 


-5.5 


4 


11.9 


IC 4663 


0.11 


3.6 


1900 


-7.7 


10 


13.4 



N V X4604, A4620 

27 Jun 

02 Jul 

03 Jul 




4620 4630 4930 4935 
Wavelength (A) 



4940 4945 4950 4955 



Figure 4. Line profiles of N V A4604, A4620 and A4944 in the 
overlap region between all three epochs. Note the absence of any 
RV shifts >5-10 km s~^ that argues against the presence of a 
post common-envelope binary companion. 



emission lines which show a root mean square radial veloc- 
ity (RV) variability of 3-4 km s~^ using cross-correlation 
techniques, i.e. no evidence of variability >5-10 km s~^. A 
main-sequence or white dwarf companion in a post common- 
envelope binary (i.e. -P < 1 day) would produce larger vari- 
ability than our detection limit and is therefore unlikely to 
be present unless the orbital inclination is very low. The 
lack of temporal coverage less than a day also makes this 
possibility difficult to rule out entirely since the median pe- 
riod of close binary central stars is 0.3-0.4 days (Miszalski et 
al. 2009a, 2011b). The small variability in the line intensity 
profile may be explained by either noise in the spectrum or 
wind turbulence (e.g. Grosdidier et al. 2000, 2001b). There 
is little historical information on any photometric variabil- 
ity of IC 4663, e.g. whether an LMC-N66-like outburst has 
occurred. The insensitivity of Shaw & Kaler (1989) to mag- 
nitudes fainter than V = 16 provides a reasonable lower 
limit to the central star magnitude in 1989. The very bright 
nebula precludes measuring central star magnitudes from 
photographic plate material. 



3.4 Variability 

The GMOS observations were taken on three separate 
epochs, namely El (June 27), E2 (July 2) and E3 (July 
3), separated by E2 - El = 5.273 d, E3 - E2 = 1.005 d 
and E3~ El ^ 6.278 d. Figure [4] depicts the narrow N V 



4 THE NEBULA 

4.1 Observations and basic properties 

The nebulae surrounding WR central stars contain valuable 
additional information to better understand the nature of 
their nuclei. Acquisition images of IC 4663 were taken during 
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Table 3. Basic nebula properties of IC 4663. 



PN G 


346.2-08.2 




RA 


17''45™28=6 




Dec. 


_44°54'17 




Nebula diameter 


16.0" X 19.5" 


this work 


Halo dimensions 


52.5" X 55.0" 


this work 


logFiHp) 


-11.44 


Shaw & Kaler (1989) 


^cxp 


30 km s-l 


Hajian et al. (2007) 


""hrv 


-78 ±3 km s-l 


this work 


tncb 


~4400 yrs 


this work 


thalo 


~30500 yrs 


this work 



our GMOS programme through the filters Ha (120 s), OIII 
(120 s) and OIIIC (90 s) whose central wavelengths/FWHMs 
are 656.0/7.2 nm, 499.0/4.5 nm and 514.0/8.8 nm, respec- 
tively. The seeing was 1.5" (FWHM) on 27 June and 1.1" 
(FWHM) on 2 and 3 July as measured from the OIIIC im- 
ages. 

Figure [5] shows the GMOS [O III] image overlaid with 
a colour-composite image made from the HST images pub- 
lished by Hajian et al. (2007). Revealed for the first time is 
a low surface brightness AGB halo (Corradi et al. 2003). A 
brightened rim lies 22" SE of the central star and is likely 
caused by ISM interaction. At ^^0.25% of the main nebula 
peak [O HI] intensity, i.e. 400 times fainter, roughly corre- 
sponding to the extent visible in the scaling selected for Fig. 
[5] the halo measures 52.5" x 55.0" across. The halo may 
even reach as far as 72" across before it becomes indistin- 
guishable from the sky background. 

The PN morphology is best captured by the aforemen- 
tioned Hajian et al. (2007) HST images that were taken on 
04 August 1998 through the i^502iV ([O HI]), F658N ([N II]) 
and F555W filters. The data were retrieved from the HLA 
and have total exposure times of 520 s, 320 s and 1020 s, re- 
spectively. A nebula diameter of 16.0" x 19.5" was measured 
at 10% of peak intensity in the F502N image. The overall 
morphology appears to be an ellipsoidal bubble surrounded 
by a second fainter shell similar to other bright PNe (e.g. 
NGC 1501). This interpretation is consistent with the high 
resolution spectrograms published by Hajian et al. (2007). 
As previously mentioned, low-ionisation filaments are also 
present (Goncjalves et al. 2001, visible as red filaments in 
Fig. [5| and were part of the selection criteria for the inclu- 
sion of this PN in our GMOS programme. 

Combining the equatorial expansion velocity of 30 km 
s~^ measured by Hajian et al. (2007), with the minor axis 
nebula diameter of 16" and our assumed distance of 3.5 
kpc (SSV08), gives an expansion age of tncb ~ 4400 years. 
For the halo we assume a 15 km s~^ expansion veloc- 
ity (Corradi et al. 2003) to find thaio ~ 30500 years and 
ihaio — tncb ~ 26000 years. We derive a heliocentric RV Uhrv 
of —78 ± 3 km s~^ from the brightest emission lines in the 
GMOS nebula spectra. This value is in agreement with —74 
km s~^ measured by Beaulieu et al. (1999). Shaw & Kaler 
(1989) measured integrated fluxes for the nebula and found 
log F{Hl3) = -11.44 (-10.93 dereddened), where F{Hl3) 
is in units of erg cm^^ s^^. Table [3] summarises the basic 
observed properties of IC 4663. 



4.2 Chemical abundances and plasma parameters 

The nebular spectrum of IC 4663 was previously analysed 
by CCMIO using low resolution (~6A FWHM) spectroscopy 
obtained with a 1.6-m telescope. Figure [6] shows our com- 
bined nebular spectrum made from averaging both of the 
2400 s GMOS exposures. We compared the measured emis- 
sion line intensities in our deep GMOS spectrum with the 
observations of CCMIO and found a good correlation with 
the line fluxes for the stronger lines, but also a tendency for 
the weak line fluxes to be overestimated by CCMIO. For in- 
stance, the He I lines at A4471 and A6678 are twice as bright 
as expected from the A5876 line, whereas all four lines in our 
spectrum are close to their theoretical ratios. Similarly, the 
fluxes of [Ar IV] A7170 and A7237 are ~3 times stronger 
than predicted from their A4740 line. Thus we shall ignore 
them in our analysis. Our [Ar IV] A4740 line is only 1/3 of 
the CCMIO value. 

Since the GMOS spectrum is deeper with greater S/N 
we produced a composite spectrum by using line inten- 
sities from our deep GMOS spectrum for the blue part 
(A s; 5900A) and CCMIO for the red part (except for [O II] 
A3728 which lies outside our GMOS spectrum). The pub- 
lished intensities of CCMIO were reddened with the redden- 
ing law of Fitzpatrick (1999) for Ry = 3.1, and the fluxes 
were scaled to the GMOS H/3 flux. We measured UV line 
fluxes from archival lUE spectra SWP33946 and LWP13707, 
observed on 20 July 1988, that were joined to the optical 
spectrum, taking into account the finite coverage of the lUE 
aperture. Although the aperture includes the central star, 
the low-dispersion lUE spectra only show nebula emission 
lines and no features from the central star. We estimate the 
stellar contribution to the lUE aperture at He II A1640 to 
be ~5% in comparison with our CMFGEN model. 

Tables 2] and [S] contain the full composite nebula spec- 
trum of IC 4663. The line fiuxes were analysed with the 
plasma diagnostics program hoppla (Acker et al. 1991; 
Koppen et al. 1991; Girard et al. 2007) which first deter- 
mines the extinction constant - using the Rv =3.1 redden- 
ing law of Fitzpatrick (1999) - as well as the electron tem- 
perature and density in a consistent way. The Hq:/H/3 ratio 
would give c = 0.58, the inclusion of the blue Balmer lines 
gives a slightly smaller value of 0.54. We adopted c — 0.51 
to be consistent with the dereddening of the stellar spec- 
trum. In either case, all Balmer lies are fit within 10% of 
their intensity from case B recombination theory. The He I 
line intensities are matched better than 10% by case B re- 
combination emissivities corrected for coUisional excitation 
(Porter, Ferland & MacAdam 2007), except for the A6678 
and A7065 lines from CCMIO. We note that the choice of ex- 
tinction constant has only very little infiuence on the derived 
elemental abundances (less than 0.04 dex). 

With the adopted extinction, we then determined the 
aperture factor for the lUE spectrum by minimizing the er- 
ror in the He II UV and optical lines AA1640, 4686, 3203, 
according to case B recombination theory. This scaling fac- 
tor is also consistent with the ratio of the nebular area and 
the lUE aperture (an ellipse measuring 10"x20"). 

Ionic abundances were derived from each of the dered- 
dened line intensities, listed in Tab. [H] using the electron 
temperatures for the low, middle, and high ionic species 
(given in Tab. [7]), and weighted with the emissivity of each 
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Figure 5. Gemini South GMOS [O III] image of IC 4663 (background) showing the newly discovered AGB halo, overlaid with an HST 
colour-composite image (red, green and blue channels were made from F658N, F502N and F555W images, respectively). The image is 
centered on the central star and measures 60 X 60 arcsec^ (~1 pc across at 3.5 kpc) with North up and East to left. 



line. Elemental abundances were obtained by summing over 
the ionic abundances and correcting for unobserved stages 
of ionisation: 

He _ He+ He++ 

IT ^ "H^^ H+ 

Carbon is derived solely from the C III] emission line (see 

below) : 

C _ C++ 
H ^ H+ 
We use the standard correction factors (e.g. AUer 1984): 



N 



N+ 0/H 



H 


H+ 0+/H+ 


O 
H " 


fO+ 0++\ He/H 
\R+ ^ H+ ) He++/H+ 


For 


neon, we use 


Ne 


Ne++ O/H 


H 


H+ 0++/H+ 



and note that the sum of Ne++ and Ne'^+ gives only a 
slightly higher abundance of 12 + log(Ne/H) — 8.17. For 
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Figure 6. GMOS nebular spectrum of IC 4663 with the main linos identified. See Tab. [5] for line intensities. 



Table 4. UV emission line intensities measured from the lUE 
spectra of IC 4633 (SWP33946 and LWP13707) scaled to H/3 = 
100. Since the lUE aperture does not cover the entire nebula, 
the fluxes were scaled as to give the best agreement between the 
intensities of the UV and optical He II lines. The measured fluxes 
in units of 10"^'^ erg en 
the flux values by 5.5. 



Line 


Flux 


Intensity 


CIV 1550 


<8 


- 


Hell 1640 


124 


558 


[NIII] 1753 


6: 


27: 


CIII] 1908 


36 


179 


[NelV] 2424 


25 


106 


Hell 2734 


10 


28 


OIII 2837* 


8: 


20: 


OIII 3047* 


11 


26 


OIII 3133* 


76 


162 


Hell 3204 


20 


41 


OIII 3299* 


18 


36 


OIII 3341* 


37 


71 



Bowen fluorescence line. 



sulphur and chlorine, we employ the ionisation correction 
factors derived from photoionisation models (Koppen, Acker 
& Stenholm 1991): 



S_ _ 
H ' 

CI 
H 



Ar 

IT 



s^ 

H^ 



1.43 + 0.196 



O^ 



0+ 



CI++ / He/H 



H+ yHe++/H^ 
Argon is treated by the correction from Barker (1980): 



Ar^ 



H^ 



+ ■ 



Ar 



,-3+ 



H+ 



+ 



Ar' 



,4+ 



H+ 



S++S'* 
S++ 



The resulting abundances are collated in Tab.[71 A few 
remarks are necessary. Our helium abundance is slightly 
higher than found by CCMIO. This is due to the inten- 
sity of the He I A5876 line being larger - by a factor of 
1.4 - in our GMOS spectrum. Due to the ionisation correc- 
tions this also results in lower oxygen and neon abundances. 
The intensity ratio of [N II] 6583/6548 (from CCMIO) of 
1.6 is substantially lower than the theoretical ratio of 2.9, 
which indicates problems with the separation of the lines 
from Hof. As to make use of all information, we use both 
lines, and we use the A5755 line from the GMOS spectrum. 
These choices influence our results on nitrogen, albeit not 
in a critical way. Taking the A6583 line only increases the 
electron temperature in the O"*" zone to 11200 K, giving a 
slightly lower nitrogen abundance of 8.22, and N/0= —0.48. 
If we took also the A5755 intensity from CCMIO - which is 
larger by a factor of 1.5 than the GMOS value - the elec- 
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Table 5. Nebular emission line intensities of IC 4663 scaled to 
H/3 = 100 as measured from our GMOS spectrum (first half) 
and CCMIO (second half). The CCMIO fluxes were obtained by 
reddening the published intensities with c = 0.58 and the Ry = 
3.1 law of Fitzpatrick (1999). 



Table 6. Ionic abundances (relative to H"*") and the emission 
lines from which they were derived. 

Ion Ionic abundance Lines used (A) 



Line 


Flux 


Intensity 


HI 3798 


4.8 


6.6 


HI 3835 


6.2 


8.5 


[Nelll] 3869 


79.4 


106.8 


Hcl 3889 


15.4 


20.6 


Hell 3924 


0.7 


0.9 


[NcIII]+HI 3967 


24.6 


32.1 


HI 3970 


11.5 


15.0 


Hcl 4026 


2.3 


3.0 


[SII] 4068 


2.8 


3.6 


Oil 4072 


1.5 


1.9 


[SII] 4076 


1.5 


1.9 


NIII 4097 


4.2 


5.2 


H<5 4102 


24.2 


30.5 


on 4119 


0.7 


0.8 


Hell 4200 


1.4 


1.7 


CII 4267 


2.1 


2.5 


H7 4340 


43.9 


51.5 


[OIII] 4363 


8.8 


10.3 


NIII 4379 


1.2 


1.3 


Hel 4388 


0.45 


0.4 


Hel 4472 


3.1 


3.5 


NIH 4511 


2.5 


2.8 


NIII 4634 


2.8 


3.0 


NIII 4641 


7.3 


7.8 


on 4649 


2.5 


2.6 


Hell 4686 


80.6 


85.3 


[ArIV]+HeI 4711 


6.9 


7.0 


[NcIV] 4714 


0.8 


0.8 


[NelV] 4724 


0.9 


0.9 


[ArlV] 4740 


2.7 


2.8 


H/3 4861 


100.0 


100.0 


Hcl 4922 


0.9 


0.9 


[OIII] 4959 


375.4 


364.0 


Hell 5412 


6.5 


5.5 


[ClIII] 5518 


1.2 


1.0 


[ClIII] 5538 


1.0 


0.8 


OIII 5592 


0.4 


0.4 


Nil 5667 


0.3 


0.2 


Nil 5680 


0.5 


0.4 


[Nil] 5754 


0.8 


0.6 


Hel 5876 


11.8 


9.0 


[Oil] 3728 


38.9 


54.6 


[OIII] 5007 


1085.0 


1035.1 


[SIII] 6312 


3.9 


2.7 


[ArV] 6435 


0.8 


0.5 


[Nil] 6548 


31.3 


20.8 


Ha 6563 


454.2 


301.5 


[Nil] 6583 


50.1 


33.5 


Hcl 6678 


5.2 


3.4 


[SII] 6716 


3.9 


2.6 


[SII] 6731 


4.9 


3.2 


[ArV] 7006 


1.9 


1.2 


Hel 7065 


2.5 


1.5 


[Arlll] 7136 


50.9 


30.8 


[ArlV] 7170 


2.1 


1.3 


[ArlV] 7237 


2.0 


1.1 


[Oil] 7325 


5.1 


3.0 


[Arlll] 7751 


11.9 


6.6 



He+ 


6.43 


X 


10-2 


5876, 


4471, 


6678, 7065, 


He++ 


6.58 


X 


10-2 


1640, 


4686, 


3203, 2734, 


C++ 


1.57 


X 


lo-'' 


1908 






N+ 


6.42 


X 


io-« 


6583, 


6548, 


5755 


N++ 


1.22 


X 


10-4 


1753 






0+ 


1.79 


X 


10-^ 


3728, 


7325 




0++ 


2.32 


X 


lo-'' 


5007, 


4959, 


4363 


Ne++ 


5.92 


X 


10-5 


3869 






Ne3+ 


8.73 


X 


10-5 


2424, 


4714, 


4724 


S+ 


5.02 


X 


10-7 


4076, 


4068, 


6717, 6731 


S++ 


4.78 


X 


10-^ 


6312 






C1++ 


7.66 


X 


10-8 


5517, 


5537 




Ar++ 


1.96 


X 


10-^ 


7135, 


7750 




Ar«+ 


3.38 


X 


10-7 


4740 






Ar4+ 


1.42 


X 


10-7 


7005, 


6434 





tron temperature increases further to 13900 K, the nitrogen 
abundance to 8.27, and N/0= —0.35, the same value as 
derived by CCMIO. Furthermore, from the just discernible 
[N III] A1753 line in the lUE spectrum, we may derive an 
ionic abundance which is substantially larger than from N+. 
Since in a high excitation nebula nitrogen is found predom- 
inantly in the form of N++, the derived nitrogen abundance 
of 12 + log(N/H) = 8.11 is a little less than that obtained 
from the [N 11] lines, but with a substantial ionisation correc- 
tion factor (30-60 in the various cases considered here) . This 
gives N/0=— 0.59. We conclude that the nitrogen abun- 
dance is certainly in the range 8.1-8.3, hence N/O = —0.6 to 
—0.3. This is just below the threshold value for a bona-fide 
Type I PN in Peimbert's classification and well below the 
(linear) value of 0.8 for a Type I in the sense of Kingsburgh 
& Barlow (1994). 

The C++ abundance derived from the optical recombi- 
nation line at 4267 A is about 30 times higher than from 
the coUisionally excited doublet C III] 1908 A. This abun- 
dance discrepancy factor is larger than the usual factor of 
2-3, but other extreme examples are known (e.g. Liu et al. 
2006) . The presence of N III and O III recombinations lines 
is in agreement with the very high level of nebular excita- 
tion, which is in concordance with the UV flux distribution 
predicted from our model atmosphere for the central star. A 
recombination line analysis of the nebula spectrum will not 
be performed here. 

Despite the high nebular excitation the lUE spectrum 
shows no detectable C IV A1550 emission. In nebulae with a 
large central hole, this line is strongly weakened, as the C*+ 
zone is remote from the central star, and hence the ioniz- 
ing radiation is geometrically strongly diluted (c.f. Koppen 
1983). The monochromatic images of IC 4663 indeed sug- 
gest a bubble-type geometry for the emission region. Tests 
with photoionisation models indicate that in a thin-shell 
type nebula the C IV lines can be suppressed sufficiently 
to agree with the lUE spectrum. As a consequence, we may 
assume that the C++ zone nearly completely fills the neb- 
ula, and that the carbon abundance derived from the C III] 
1908 lines may be close to the true value. Since the neb- 
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Table 7. Plasma diagnostics and the chemical composition of 
IC 4663, given in the usual notation 12+log[n(X)/n(H)], from this 
work and from CCMIO. For comparison the average abundances 
for Type I and non-Type I PNe are taken from Kingsburgh & 
Barlow (1994) and Solar abundances from Asplund et al. (2009). 





IC 4663 


Type I 


non- 


Sun 




this work 


CCMIO 




Type I 




c 


0.51 


0.58 











ne [cm"-^] 


1000 


3810 


— 


— 


— 


Te(0+) [K] 


10300 


12800 


— 


— 


— 


Te(0++) [K] 


11500 


11500 


— 


— 


— 


Te(He++) [K] 


11400 


— 


— 


— 


— 


He 


11.11 


11.07 


11.11 


11.05 


10.93 


C 


8.20 


— 


8.43 


8.81 


8.43 


N 


8.26 


8.49 


8.72 


8.14 


7.83 


O 


8.70 


8.84 


8.65 


8.69 


8.69 


Ne 


8.11 


8.25 


8.09 


8.10 


7.93 


S 


7.02 


7.03 


6.91 


6.91 


7.12 


CI 


5.50 


— 


— 


— 


5.50 


Ar 


6.43 


6.87 


6.42 


6.38 


6.40 


log(C/0) 


-0.50 





-0.12 


+0.12 


-0.26 


log(N/0) 


-0.45 


-0.35 


+0.07 


-0.55 


-0.86 


log(Ne/0) 


-0.59 


-0.59 


-0.56 


-0.59 


-0.76 


log(S/0) 


-1.68 


-1.81 


-1.74 


-1.78 


-1.57 


log(Cl/0) 


-3.20 


— 


— 


— 


-3.19 


log(Ar/0) 


-2.27 


-1.97 


-2.23 


-2.31 


-2.29 



ula has a rather elhptical shape and also shows indications 
for a double shell, we refrain from constructing a spherical 
photoionisation model. 

The derived abundances show that our He/H is slightly 
higher than in CCMIO and closer to a Type I composition 
(Peimbert & Torres-Peimbert 1983; Kingsburgh & Barlow 
1994), because the He I 5876 intensity from the CMOS spec- 
trum is higher. We note that our value is also consistent with 
the other He I lines. However, the N/0 ratio is not as high 
as to classify it as a genuine Type I nebula. The nearly So- 
lar oxygen abundance is also similar to the average PN, as 
are the Ne/O and S/0 ratios. The ratios of the a-elements 
are close to the average values found in H H regions of spi- 
ral and irregular galaxies (Henry & Worthey 1999), namely 
Ne/O = -0.67 (slightly lower than in PN and in IC 4663), 
S/O = -1.55 (slightly higher than in IC 4663), and Ar/0 
= —2.25. Since these elements are not expected to be af- 
fected by nucleosynthesis in the progenitor star, their ratios 
reflect the interstellar medium values which are determined 
by nucleosynthesis in massive stars. 

The carbon to oxygen ratio is lower than in both PN 
types and the Sun, but in light of both the strong sensitivity 
of the UV lines to electron temperature and the uncertainty 
of the ionisation correction factor for carbon, this difference 
of 0.3 dex cannot be taken as a solid proof for a carbon 
underabundance. If placed in Fig. 7 of Henry & Worthey 
(1999), IC 4663 would be located near the lower envelope 
of scattered points, but would still be consistent with the 
average relation. Finally, the chlorine to oxygen ratio is So- 
lar. The abundance pattern of IC 4663 is close to the Sun, 
except for enhancements in helium, nitrogen and neon. 
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Figure 7. Hydrogen-deficient PNe in the distance-independent 
log Sv and log Sh/} plane (see Corny & Tylenda 2000). 



4.3 Evolutionary status 

The evolutionary state of the PN IC 4663 may also be com- 
pared against other PNe using the distance-independent pa- 
rameters log Sv and log Snp (Gorny & Tylenda 2000). The 
parameters are defined as the H/3 surface brightness and an 
analagous quantity replacing Fnp with Fv, an extinction- 
corrected stellar flux in the 1/-band, respectively. Both Sv 
and Sh/s decrease as a PN expands and evolves towards the 
WD state. Figure[7]portrays the position of IC 4663 amongst 
other PNe with H-deficient central stars (log Sv = —6.54 
and log Shi3 = -2.71) and PB 8 (log Sv = -4.11 and log 
Shp = —1.48, based on the updated parameters in Todt et 
al. 2010a). The location of IC 4663 is coincident with PNe 
that have a [WO] central star (i.e. [WC4] or earlier) and 
PB 8 with those that have a [WC] central star (i.e. [WC7] 
or later). 



5 DISCUSSION 

5.1 A bona-fide PN 

Table [8] delineates three illustrative distance scenarios for 
IC 4663. At an assumed distance of 3.5 kpc (SSV08), the 
central star has My — +3.1 mag. This is over 6 magnitudes 
fainter than Afy = —3.0 mag, the typical luminosity for a 
massive WN3 star (Crowther et al. 1995b; Hamann et al. 
2006). The other two scenarios demonstrate that the sub- 
luminous nature of the central star is robust against any 
reasonable uncertainties that may be associated with the 
distance of SSV08. Even if IC 4663 were a member of the 
Galactic Bulge at d = 8.0 kpc, then the central star is still 
at least 4 magnitudes too faint. IC 4663 cannot however be 
a member of the Bulge because its angular diameter and H/3 
flux are anomalously large compared to most Bulge PNe (see 
Fig. 7 of Corny et al. 2009). Indeed, an implausible distance 
of 58 kpc would be required to reproduce AIv ~ —3.0 mag, 
therefore for any reasonable distance estimate the central 
star must be sub-luminous and belong to the PN. At 3.5 
kpc IC 4663 would lie 500 pc below the Galactic Plane, i.e. 
between the bias-corrected scale heights of the thin (300-375 
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Table 8. Absolute visual magnitudes and vertical displacement 
from Galactic Plane for the central star of IC 4663 at a range of 
distances {V = 16.90 mag and E{B -V) = 0.35 mag, Sect. E^J. 



d 
(kpc) 



A/v 



Comment 



(kpc) 



3.5 +3.1 -0.50 

8.0 +1.3 -1.14 

58.0 -3.0 -8.27 



d from SSV08 

d of Galactic Bulge PN 

Mv of massive WN3 star 



pc) and thick (900-940 pc) disks (Siegel et al. 2002; Juric 
et al. 2008). Massive WR stars that are not runaway stars 
would be restricted to the young disk population which has a 
scale height of only ~ 48 ± 3 pc (Bonatto et al. 2006). This 
intermediate height might on its own suggest a thick disk 
progenitor, however the solar oxygen and argon abundances, 
which reflect the metallicity of the gas which formed the 
progenitor star, suggests IC 4663 is a thin disk object. 

The presence of an AGB halo is a telltale signature of 
a PN (Corradi et al. 2003). The ellipsoidal morphology of 
the main nebula, its relatively low expansion velocity and 
density are all typical of PNe. Its evolutionary status over- 
laps with other PNe in Fig. [7] and exhibits a highly ionised 
emission line spectrum with a typical chemical abundance 
pattern for a PN. 



5.2 A second post- AGB H-deficient and He-rich 
evolutionary sequence 

Table |9] contrasts the He-dominated composition of the 
[WN3] central star of IC 4663 against those of other 
H-deficient post-AGB stars|j Clearly IC 4663 does not 
fit the carbon- and oxygen-rich evolutionary sequence 
[WCL]^[WCE]->PG1159 (Werner & Herwig 2006). Rather, 
it fits amongst other rare He-dominated objects includ- 
ing RCB stars, 0(He) stars, extreme helium B stars and 
helium-rich subdwarf O stars for which there is no ac- 
cepted explanation for their origin (Werner & Herwig 2006). 
The born-again evolutionary models simply cannot repro- 
duce the helium-dominated atmospheres (e.g. Herwig 2001). 
Some proposed explanations include a double-degenerate 
merger (Saio & Jeffery 2002), a low-mass star leaving the 
AGB shortly after the first thermal pulse (Miller Bertolami 
& Althaus 2006) or a diffusion induced nova (Miller Berto- 
lami et al. 2011). 

If we speculate that IC 4663 has a companion, then 
perhaps a close binary evolution channel may be responsi- 
ble for the helium-rich composition. Yungelson et al. (2008) 
proposed such a scheme to explain helium-rich AM Canum 
Venaticorum (AM CVn) stars (Warner 1995; Solheim 2010). 
IC 4663 may be the ejected common-envelope at one of 
these stages making it a possible AM CVn progenitor. If 
IC 4663 were a pre-CV, then it could potentially evolve into 
a helium-rich nova after additional angular momentum loss 
(e.g. V445 Puppis, Woudt et al. 2009). Another possibility 
is that a common- envelope merger already took place and 



* Some caution should be taken when comparing compositions 
derived from different analysis techniques (see e.g. Werner & Her- 
wig 2006). 



this would be consistent with the apparent lack of evidence 
for a close binary (Sect. l3^ . More generally this would help 
explain the apparent lack of [WR] central stars in close bina- 
ries, despite many nebulae around [WR] central stars show- 
ing post-CE morphologies (Miszalski et al. 2009b). 

The discovery of the [WN3] central star of IC 4663 
has several implications. Joining the 0(He) stars (Mendez 
1991; Rauch et al. 1994, 1996, 1998, 2008) it is now the 
fifth hot (Tcft > 100 kK), He- rich and H-deficient post- 
AGB star known. A second post-AGB H-deficient sequence 
[WN]^-0(He) now seems to be a logical conclusion based on 
the strong match between their compositions. This has been 
hypothesised before (e.g. Werner 2012), but only now does 
it seem to be real with the [WN3] discovery in IC 4663. The 
RCB stars may also play an uncertain role in this sequence 
as RCB^'[WN]^-0(He). The nebula evolutionary status of 
IC 4663 and the two PNe with 0(He) central stars (Sect. 
I4.3|l . being coincident with [WO] and PG1159 stars, respec- 
tively, supports the existence of a second sequence as their 
evolutionary status mirrors the [WC]— >-PG1159 sequence. 
Once the stellar wind of the [WN3] central star of IC 4663 
disappears it will likely become an 0(He) central star and 
later a DO WD. 

In summary, we propose the [WN]— >-0(He) sequence 
is a logical conclusion from the strong match between the 
surface compositions of the [WN3] central star of IC 4663 
and the 0(He) stars (Tab. [§}. The evolutionary status of 
their nebulae (Sect. I4.3|l suggests it forms a second paral- 
lel H-deficient evolutionary sequence to the existing carbon- 
rich [WC]^-PG1159 sequence. The relatively large amount 
of residual hydrogen in PB 8 suggests it may belong to a 
separate [WN] sequence to IC 4663, perhaps refiecting the 
dichotomy between H-deficient early-type WN stars and H- 
rich late-type WN stars (e.g. Hamann et al. 2006). 



6 CONCLUSIONS 

We have discovered the first bona-fide [WN] central star of 
a PN in IC 4663. Our main conclusions are as follows: 

• The Gemini South GMOS spectrum of the central star 
is dominated by broad He II and N V emission lines, that 
together with the notable absence of He I and C IV, corre- 
sponds to a [WN3] subtype. 

• A CMFGEN NLTE model atmosphere of the central star 
was built and reasonably reproduces the GMOS spectrum. 
The central star has a very hot (T* = 140 kK) and fast 
("Uoo = 1900 km s"'^) stellar wind. The atmosphere consists 
of helium (95%), hydrogen (< 2%), nitrogen (0.8%), neon 
(0.2%) and oxygen (0.05%) by mass. Such a composition 
cannot be explained by extant scenarios that seek to explain 
the H-deficiency of post-AGB stars. 

• The helium-rich atmospheric composition strongly links 
[WN] central stars to the 0(He) central stars (Rauch et 
al. 1998). The data support a second post-AGB H-deficient 
evolutionary sequence [WN]— ^O(He) parallel to the estab- 
lished carbon-rich [WCL]^[WCE]^PG1159 sequence. Con- 
nections to other He-rich/H-deficient stars may include RCB 
stars, and if there were a binary companion to the [WN3] 
star, AM CVn stars and He-rich novae after further angu- 
lar momentum loss. The existence of this second sequence 
strongly suggests another mechanism exists for producing 
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Table 9. Stellar compositions expressed as mass fractions (%) of IC 4663 and other objects with similar compositions. 



Object 


Type 


Tcff 
(kK) 


logs 

(cgs) 


M 


H 


He 


C 


N 


O 


Ne 


Reference 


IC 4663 


[WN3] 


140 


6.1t 


1.8 X 10-*^ 


<2 


95 


<0.1; 


0.8 


0.05 


0.2 


this work 


K 1-27 


0(He) 


105 


6.5 


7.9 X 10-1" 


<5 


98 


<1.5 


1.7 


- 


- 


Ranch et al. (1998) 


LoTr 4 


0(He) 


120 


5.5 


2.0 X 10-» 


11 


89 


<1 


0.3 


<3 


- 


Ranch et al. (1998) 


PB 8 


[WN6/WC7] 


52 


- 


8.5 X 10-* 


40 


55 


1.3 


2 


1.3 


- 


Todt et al. (2010a) 


RYSgr 


RGB 


7.25 


0.7 


- 


6x10-'' 


98 


0.7 


0.3 


0.09 


- 


Asplund et al. (2000) 



t Adopted value assuming M = 0.6 Mq. 

: Model constrained (no carbon lines observed). 



H-deficient post-AGB stars besides the various horn-again 
scenarios. 

• For aU reasonable distances the luminosity of the central 
star is 4-6 mag fainter than massive WN3 stars (Crowther 
et al. 1995b). At our assumed distance of 3.5 kpc (SSV08), 
the derived luminosity is log L/Lq = 3.6, the stellar radius 
R* = 0.11 Rq and the mass loss rate M = 1.8 x 10"* Mq 
yT~^. The transformed radius, Rt, of 13.4 Rq for IC 4663 
is comparable to massive WN3 stars, indicating a similar 
stellar wind density to its massive counterparts. This favours 
our morphological [WN3] classification of the central star of 
IC 4663, rather than alternative weak wind (Of or Of-WR) 
classifications (Mendez et al. 1990; Mendez 1991). 

• The surrounding nebula is a bona-fide PN. The ellip- 
tical main nebula has a relatively low expansion velocity of 
30 km s^^ and is surrounded by a newly discovered AGB 
halo. The evolutionary status of the highly ionised nebula 
corresponds to that of similar nebulae around [WO] type 
central stars and supports the [WN]— 7-0(IIe) evolutionary 
sequence. 

• Chemical abundances of the nebula were derived from 
a combination of lUE, CMOS and CCMIO emission line 
intensities. The abundance pattern is approximately Solar 
with slight enhancements for helium, nitrogen and neon. The 
carbon to oxygen ratio is lower than in the Sun, however the 
data are insufficient to prove a carbon underabundance. 

• We found no evidence for RV variability from the cen- 
tral star >5-10 km s"^ in three epochs. Further monitoring 
of the central star should be performed to conclusively rule 
out the presence of a companion to the [WN3] central star. 
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